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THEORY OF THE LANDING IMPACT OF SEAPLANES.* 
By Vvilhelm Pabst. 

The present investigation is an endeavor to express the 
jolting stresses, designated as landing impacts, undergone by 
seaplanes in landing and t eking off from rough water, as func- 
tions of specific factors, in order to enable the 6 valuation of 
empirically obtained results and thus acquire theoretical data 
for the construction of seaplane floats and hulls. A physical 
explanation of the landing impact, on which the following math- 
ematical investigation is based, was given by F. Secwald,** 
Director of the Aerodynamic Section of the D.V.Li, who first 
suggested the present investigation. 

General Considerations on Landing Impacts 

1. An ideal landing on smooth water resembles a take-off 
from smooth water, except that it takes place in the reversed 
order. The bottom of the seaplane is gradually submerged when 
the step approaches the smooth water surface t anient i ally . The 
increasing drag gradually reduces the speed of the seaplane, 

simultaneously replacing the lift of the wings by the buoyancy of 

♦"Theorie des Landestosses von Seef lugzeugen. 11 From Zeitschrift 
fur Flugteclinik und Motorluf tschif f ahrt , May 14, 1930, pp. 217- 
226. 

♦♦Discussion of F. Z. Diemar f S lecture on "Flugboot und Seegung," 
1927 Yearbook of the W»Gr.L. 
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the water. The braking resistance consists of waVe-pr oduc ing 
and frictional resistance. Another resistance is due to the 
fact that , when a moving body is plunged into water, a certain 
mass of the water is accelerated, corresponding to the flow 
about the body. When the bottom plunges very quickly into the 
water, as may occur in alighting on rough water and meeting a 
wave head-on, the last-named resistance becomes so great that 
the first two can be neglected. Froude's and Reynold^-' laws 
of similarity may therefore be disregarded as applied to the 
impacts actually considered in this connection. The following 
considerations prove that Newton's law of similarity does not 
fully apply to all cases. 

2. As shown by F. Seewald, the magnitude of the impact is 
affected by the elasticity of the airplane, as well as by the 
usual factors of current methods of calculation, such as the 
roughness of the water, the weight of the seaplane, the landing 
speed and the size and shape of the float bottom. With a stiff 
float bottom o.nd assumed compressibility of the water, the 
acceleration of a finite mass of water would take place in an 
infinitely short time, thus producing, according to the momentum 
theorem, infinitely great forces. The fact that the forces do 
not become infinitely great is chiefly attributable to the in- 
fluence of the elasticity of the seaplane, since the compressi- 
bility of water is negligible in comparison. 

This consideration is confirmed experimentally by the fact 
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that the impact acceleration decreases from the float bottom up. 
This fact, which is taken care of, in the loading conditions of 
the D.V.L., by a higher load factor of the float, can be ex- 
plained only by the effect of elasticity. The process therefore 
involves elastic forces in addition to inertia forces* Hence 
Cauch^s law of similarity is applicable in this case. I T j states 
that two phenomena are dynamically similar only when Caucht ! s 



number C = v yp/E has the same value in both cases. The fact 
that Cauchy's number contains the velocity v, leads directly 
to the surmise that the landing impact does not, or not in all 
cases, depend on the square of the landing speed, as would fol- 
low from Newton 1 s law of similarity without modification. It 
should also be investigated as to whether, considering the ef- 
fect of elasticity and the short duration of the impact, the 
latter is also affected by the damping action of the material* 

3. The area of the float bottom coming simultaneously into 
contact with the surface of the water has a decisive effect on 
the magnitude of the landing impact. The size of this area de- 
pends on the region in which the water surface a&d float bo t torn 
are parallel when coming into contact with each other. These 
conditions are greatly affected by the roughness of the water 
and the manner of alighting. Certain simplifications are re- 
quired by the great variety of the seaways, which change in wave 
shape ana length according to the force ana duration of the wind, 
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the length of the unobstructed wind path and the depth of the 
sea. In mpst cases a seaway consists of several superposed sea- 
ways. Such diagrammatic representations, as axe successfully 
used for strength calculation in shipbuilding practice, hold good 
only in a very few cases but enable a mathematical estimation 
of the various factors, a comparison between the different mod- 
els and, an improvement in the load conditions, when supplement- 
ed by experimentally determined values. Hence, the following 
calculations apply to the so-called ''established* 1 seaway - a sea- 
way which, according to numerous observations,* is constant 
when there is a long enough unobstructed wind path, sufficient 
depth of water, and a steady wind of constant direction and 
force. The adopted wave shape is trochoidal. 

Similar simplifications and assumptions should be made on 
the manner of landing. It would then be possible to work out 
landing cases for specific seaways and to calculate the landing 
impact by the method set forth below, using the bottom contact 
areas obtained from the drawing. The genera.! conditions and the 
extent to which the pilot can be expected to avoid very rough 
landings are chiefly matters of opinion and depend moreover on 
the aerodynamic properties of the aircraft. Also, the seaworth- 
iness requirements vary in the different cases according to the 
structural problems to be solved. Therefore, the methods of 
calculation presented below- should be applied with discrimination 

in the different cases . , — 

♦Zimmerman, Von Koopen and Laos. Sec also Johow Foerster, "Hilfs- 
buch f.d. Schiffbau, » 5 th edition, Vol. I, pp. 41S and 420. 
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Among the different possibilities of alighting on rough 
water the maximum landing impact which a seaplane should be 
capable of withstanding is produced in the following case: The 
seaplane strikes the water flatly with the straight portion of 
the bottom. This may happen when, in attempting to land with 
the tail down, the seaplane meets an oncoming wave, or in start- 
ing, takes off prematurely and falls back on the water. The 
alighting moment is shown diagrammat ically in Figure 1. Seaway 
2 with a wave length of 11 rn (36 ft.) and a wave depth of 1 m 
(3.23 ft.) is roughly represented in the figure. The wave is 
represented as usual by a trochoid. The wind speed correspond- 
ing to the seaway is approximately 3 m/s (9.8 ft. /sec). Let 
c a be the seaplane speed with respect to the water. It consists 
of the corresponding components of the speed above ground com- 
bined with the velocity of the water. In general c can be 

a 

replaced by the corresponding component of the landing speed nor- 
mal to the keel at the maximum angle of attack. Inasmuch as the 
seaplane is usually brought down against the wind, the reduction 
in the landing speed caused by the head wind is balanced by the 
opposite motion of the water. Of course, the actual values may 
be used in any particular case. 

The Accelerated Water Mass 

4. As already mentioned above, all the other forces are 
assumed to be small in comparison with that of the impact of the 
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float bottom on the water. According to Laftfe* the following 
formulas ere then obtained under the sole action of the impulsive 
pressures of the float bottom in a two-dimensional solution, 
when the motion is started from the position of rest 

u - - i a_g 
p a x 

pay 

in which 8f is the impulsive pressure. 

By means of the equation of continuity, we then obtain 

d & -i- q Sti - 0 
By 2 

provided p is constant which, on account of the slight compress- 
ibility of water, seems admissible even for very large impact 
forces. We shall now consider the marginal conditions for a 
plate of infinite length and width b lying on the water. 

If the impulsive pressure of the above equations is replaced 
by $ = p $, in which $ is the velocity potential, and it is 
considered that no impulsive pressures are exerted on the open 
water surface, i,e« , £> = 0 and $ = 0, the same marginal con- 
ditions prevail as for the plate of infinite length on an infi- 
nite liquid surface. The following formula is then obtained for 
the plate velocity at the end of the impact period, provided the 

one-sidedness of the process is taken into consideration. 
*Lamb-Friedl , "Lehrbuch der Hydrodynamik, " paragraph 12. 
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8 J 

u = ~ 

p n b A a 

This applies to a plate element of width b and length Aa 
when J as p At . Hence the mass of water to be accelerated is 

M w = p g b 2 A a (1) 

The distribution of the impact pressures over b is elliptical. 
It is assumed that the plate is absolutely rigid. In practice 
the flow and the pressure distribution are subject to variation, 

5, The assumption of an infinite length does not apply to 
the actual float bottom. In fact, the length of the bottom 
portion which strikes the water is of the same order of magni- 
tude as the width. Since, on the assumption of an infinite 
plate length, the bottom width goes into the second power, while 
it has a smaller power in the case of a finite bottom length, 
the latter must be taken into consideration, in order to avoid 
wrong conclusions regarding the effect of the width of the hull 
on the impact. I n this case the bottom portion concerned can 
also be considered as a plate in an infinite liquid, the one- 
sidedness of the process being taken into consideration. The 
water mass accelerated by such plates of a finite chine ratio 
was determined experimentally by means of small vibrations. Whe 
a body vibrates in a nonviscous, incompressible, infinite fluid 
at rest, the mass of the body is increased by the flow which de- 
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velops during the motion ( Stokes* ). Green** ) . In an ideal non- 
viscous fluid this flow is a potential flow. In viscous fluids 
the potential flow can he maintained with a good approximation, 
provided the body makes very short and quick vibrations.*** 
Under the above assumptions of disregarded friction and wave 
formation, this fact permits of easily determining the acceler- 
ated water mass as already suggested by Fottinger for other pur- 
poses but, so far as I know, never put into practice. 

Figure 2 shows the test installation. Plate 1, stiffened 
by a longitudinal rib, is secured to a duralumin tube 2, which 
connected with two steel springs 3 and 4, can vibrate along 
its longitudinal axis. This system, which is capable of vibrat- 
ing, is deflected approximately 0.2 mm (0.008 in.) and then sud- 
denly released by severing a wire. The resulting damped vibra- 
tion was plotted by means of a scratch recording device 5,**** 
directly and without lever transmission, with a diamond on a 
glass plate moved laterally by the electric motor 6. The re- 
sulting diagram was estimated under a microscope using the simul- 
taneously recorded time marks. This estimation, made on the as- 
sumption of proportional damping, showed that the influence of 
the damping on the period of vibration Was negl igibly small. 

*Stokes, "On Some Cases of Fluid Motion." Camb. Trans., 8, 
1843, Math, and Phys. Papers I, p. 17. _ 

**Green, "Researches on the Vibration of Pendulums m Fluid 
Media." Trans. R. S. Edin. , 1883, Math. Papers, p. 315. 

***Fottinger, Jahrbuch d. Schif fbautechn. Gesellschaf t , 1924. 
****Pabst, W. , "Aufzeichnungen schneller Schwingungen nach dem 
Ritzverfahren,* Zeitschrift des Vereines deutscher Ingenieure, 
1929. No. 46. ' 
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Figure 3 is a microphotograph of such a vibration and of the cor- 
responding time marks. Before the tests, the spring constant 
was determined by loading the device with known weights and re- 
cording the resulting deflection (Table I). 



TABLE I. Determination of Spring Constant 



No. 


Load 
kg 


Distance from 
base line 
l/lOO mm 


kg/ cm 


Mean K 
value, 
kg/ cm 


1 


2.5 


14.3 


175.0 






2 


3.0 


17.0 


176.5 






3 


3.5 


20.0 


175.0 






4 


4.0 


23.0 


174.0 






5 


4.5 


25.5 


176.5 






6 


5.0 


28.0 


178.5 




176.0 


7 


5.5 


31.0 


177.5 






8 


6.0 


34.0 


176.2 






9 


6.5 


37.0 


175.9 






10 


7.0 


40.0 


175.0 






11 


7.5 


42.5 


176.5 















The mass of the instrument was then determined by causing it to 



vibrate in air. A comparison of the mass determined by vibra- 
tion with that obtained by weighing showed that the steel 
springs participated in the vibrating mass of the device to an 

extent of 35.8$ of their total mass. The vibration of the plates 
(the dimensions and weights of which are given in Table II) 

against water was then tested by placing the device over the 
water-filled tank shown in the background of Figure 2. The water 
surface w a s approximately 45 cm (18 in.) above the plate and did 
not seem to be affected by its vibration. At the point of im- 
mersion of the tube, a concentrically progressing undulatory 
motion of very small amplitude was observed. It was merely due 
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to friction and capillary action and probably did not affect 
the vibration of the plate. The longer of the four plates test- 
ed were stiffened in the plane of symmetry of the flow in order 
to avoid natural vibrations of the plate. The plate edges were 
rounded off elliptically . The results are given in Table III. 



TABLE 


II. Dimensions 


and 


Weights of 


Plates 


Plate 


a 


b 


a:b 


G 


Ho. 


mm 


mm 




g 


1 


100 


100 


1 


82 


2 


200 


100 


2 


247 


3 


500 


100 


3 


343 


4 


400 


100 


4 


460 



All moved parts without plates and springs G = 170 g 
Weight of springs G = 420 g 



TABLE III. Results 



Plate 
Ho. 


T 
s 


Total 
mass 

M 


Mass of 
device 
M A 


Mass of 
wat er 








** 
cm 


cm 


S^ 

cm 


— s 2 
cm 



1 against air 0.0096 - - - - 

1 " water 0.01391 0.863 0.410 0.453 0.453 

2 » ■ 0.01990 1.752 0.578 1.174 0.146 

3 11 11 0.02438 2.632 0.677 1.955 0.073 

4 " ■ 0.02340 3. 590 0.796 2.79 4 0.050 

The four measured values provide seven points of the diagram in 
all, since the moved mass can be easily calculated for the recip- 
rocal edge ratio. The results are plotted in Figure 4. 

As was to be expected, the curve runs, for large a/b val- 

TT 

ues, parallel to the line P q ^ a for a plate portion of the 
length a of the infinitely long plate. Thus, for the finite 
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plate length, a constant value p y^r 15 3 , must be deducted for 
the edge effect when ^ 5 1 , so that 



For very small values the curve approaches the parabola 

p — a 2 b. Approximately within the range Jr » 1 - 3, which con- 
o b 

cerns the present problem, a material reduction of the acceler- 
ated v/ater mass is achieved by taking the finite edge ratio into 
consideration. 



6, The whole seaplane-float system is based on Figure 5. 
Let the mass U x of the seaplane bu concentrated in one point 
and a spring exerting a force P = kf, assumed to have no mass, 
be fitted between the float bottom also, assumed to be without 
mass. A certain, water mass M a is accelerated when the flat 
float bottom strikes the water. 




(2) 



Impact of a Flat-Bottomed Seaplane 



Mass of seaplane 



Mass of accelerated water 



M 2 



Force of spring 



F 



kf 



Deflection of spring 



f 



Length of spring 



L 



Effective weight after 
deducting wing force 
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Then 



For 



d 2 e ^ 

M , — 4- = kf - v G I 
1 d t 2 



M a § = - kf 

d t 2 

X, - Xn = L - f. 



d 2 x, d 2 x a d 2 f 
L + 1 = - , 

d t* 



df d V 



d 2 f 



/I 1 N 
+ kf f — * — > - wg = 0. 

d XT, v Mj l£ s y 



The solution is 
where 

for 
and 

whence 



f = A sin a) t + B cos u) t + 



0)' 



a) 2 



. M„ + M. k 

f = 0 
a) t - 0 



B = - 



The constant A results from the following consideration. 
Since x x represents the path of the C.G. of the seaplane, 

2*1 = c- for cot = 0, 
dt "* 

x 2 is the path of the water mass. For oo t = 0 

*3l = o, 

dt 
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Hence 

d f d x 2 d x 1 

d~~t = dfl d "fc = " ° a? 

v/hence 

A CO 

The solution therefore reads 

t « sin co t — cos a> t + 

00 0) 2 O) 2 

= A sin ((o t - (3) + -ii, 



where 
said 



oo v \ to y a ^ 



3 - arc tan 

C 



for the maximum impact force 



P = * fmax = A + ^ =7 (i> g u) 2 + c a 2 k y + v g |i (4) 

where 



This formula is quite general and therefore applicable to 
landplanes as well, provided the mass of the wheel is disregard- 
ed. In this case M 2 =? oo and the expression becomes jjl = M x . 
When c a = 0 and V = 1, then P » S 3, which corresponds to 
Poncelet ! s theorem of mechanics. For a seaplane M 2 is usually 
only a fraction of the mass. 

The water mass obtained by approximate calculation for the 
HE 5 Heinkel monoplane corresponds to the straight portion of 
the flat bottom and is 20$ of the mass of the seaplane. Thus 
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tl 



[\e expression 



in the formula is small in comparison 



M 



+ M. 



with the expression c a 2 k u, especially when it is recalled 
that the spring constant k of seaplanes is very large and that 
the lifting forces have not, in most cases, become zero, so that 
V «* It We can therefore neglect the expression v g |i and ob- 
tain, for the landing impact, 



On the preliminary assumption that the impact is propor- 
tional to \ 2 = (i/Lj 2 Where I in an arbitrary length of one 
and L the corresponding length of another larger and geomet- 

p 7 X 3 

rically similar float, the stress a = s cgr = 1 is there- 

W A, 

fore independent of the increase in length, a condition which 
should be required. For the same material -and modulus of elas- 
ticity E, the spring constant k is proportional to \ . 
Hence, since M is proportional to A 3 , B rArjX ** X 2 > as above 
temporarily assumed. 

The impact is therefore proportional to the dimension of a 
surface and the impact load P/F is independent of the size of 
the seaplane, provided geometrical similarity is assumed. More- 
over, the elastic impact represents a vibration of the system, 
its mass consisting of the bodies which strike each other. The 
initial conditions of the vibration are determined by the rela- 



P 



max 




where 




(5) 
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tive motion of the two bodies at the beginning of the impact. 
The eccentric impact for the considered approximation is 
obtained in the usual way when the reduced mass is substituted 
for the actual mass of the seaplane 

i + r 2 

where i is the inertia radius and 

r the distance of the percussion force from the.C.G. 

7. The above approximation method is not quite satisfac- 
tory, since the seaplane itself is assumed to be rigid, while 
only the bottom is considered elastic. The same load factor for 
all parts is therefore obtained by calculation. As a matter of 
fact, the elasticity is distributed over the whole mass system, 
so that the seaplane portions nearest to the point of applica- 
tion of the impact forces, such as the floats, must withstand 
greater impact forces than more distant parts, such as the fuse- 
lage, wings and engines. 

Yet an accurate calculation seems impossible. Therefore 
we must endeavor to divide the whole system into separate masses 
connected with each other and with the water mass by elastic 
members having no mass. Analytically speaking, the separation 
of the seaplane into two masses leads, in general, to difficult 
calculations, but supplies analytically simple results which 
make it possible to answer a number of questions. The calcula- 
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tion is based on the diagram in Figure 6. Assume the fuselage, 
with the engine, wing and tail surfaces, to be rigid and have a 
mass M 1 and a moment of inertia 0 1 and represent it by the 
line AB. The float, likewise assumed to be rigid, is given a 
mass M s and represented by the line CD. The float and fuse- 
lage are connected by two elastic members which are assumed to 
have no mass. The spring constants of these members are k v and 
kj|« Between the float and the water mass there is also insert- 
ed an elastic member which represents the float bottom and which 
has a spring constant k 2 . The other notations are given in 
Figure 4. 

The manner of calculation is derived from Lagrange's equa- 
tion. 



d 



a e 

d t 



9 E 

a Qi 



in which 



E is the energy of the whole system, 



the coordinate of the respective C.G 




the first derivative after 



time 



t 



the external force. 



Hence 



E 




"a f d *a > 8 
2 Vd t y 
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Then, for the mass M x , we have in the X direction 

d 9 E 



d t 1 d t 2 ' a x : 

Moreover, the elastic forces are introduced as an external force, 
so that 

Fi = ky [(x x + a (p x ) - (x 2 + L 1 + a cp 2 )] + 

+ k h [x, - b <p a ) - (x 3 + Lj - b cp 2 )] . 

Then „ 

M x x + (k v + k h ) f 12 + (k v a - k h b) cp = 0, 

when 

x i ~ X 2 = ^i ~ *13 
cp, - cp 2 = cp 

Similarly, for q = cf> x , 

6 L cp + (k v a - k h b) f 12 + (ky a 2 + k h b 2 ) cp = 0 

the 

and for/ other masses M 2 G 2 with their coordinates x 2 and 
cp 2 and M a with x 3 . The following expressions are then ob- 
tained for the present problem 



d 2 x 

M g j-p- -(k p + k h )f 12 -(ku a - k h b)cp 12 + k, f 23 = 0 
d 2 cp 

% cf~F " (ku a " kilb)f ^ ~ (ky a2 + k h t)2 ) f * ) i2 + k 2 r f 23 = 0 

M 3 |^|f - k 2 f 23 = 0. 



N.A.C.A. Technical Memorandum No. 580 



18 



After a double differentiation of the formulas 
X x - x 2 = L ~ f xs and x 2 + rcp 2 - x 3 = L 2 - f S3 and c? - <p 2 = cp 

These equations can be introduced into the system of equations 
for the eccentric impact of the two-mass system for any arbitra- 
ry distribution of elasticity between the two masses. 



a t 



d 2 cp 



1 2 



d f 



d 3 f. 



23 



d t : 



Mi M 2 / 

+ (kp a - k h b) (— + cp i2 - — ^ 3 



- (ky a - k h b) (~ + — ) f ls + 



2 3 



i /! r2 1 \ * 

/feu + kh + a - k h b ^ 



/ 12 
2 



/kp a - b k u a + k* n b \ 
f H — ■ + „ £J r i cp 



0- 



7) 



/ 12 ^ 



As shown by the preliminary formulas, the problem leads to 
vibrations with force combinations similar to the torsional vi- 
brations of shafts. The graphical methods of calculation used 
for torsional vibrations of multipass systems may possibly be 
advantageously applied to the present case. It should be taken 
into consideration, however, that the problem differs slightly 
in the present case* From the analytical viewpoint, the solu- 
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tion of this system of equations requires a very complicated 
calculation. 

8. For the centric impact (step impact), r = 0- More- 
over, the equations with 

kv> = kv = — and a - b = — 

are simplified whence there is obtained for the impact force 

M, - p, = f 13 
and M 2 - P 2 = kg f s3 

where f 12 and f 23 are derived from the system 

k. 



d 2 f. 



i_2_ 



d V 



+ k , — + — > f. a - f 23 = 0 
1 M s -' 12 M 2 23 



d 2 f 



(8) 



d t 



3 2 v M 2 M 3 / 23 M 2 12 



or 



d 2 f. 



i § 



d t 
d 2 f . 



+ a i 2 f 12 - a s 2 f 



23 



0 



33 



d t" 



+ a 3 2 f 23 - a 4 2 f ia = 0 



This system can be easily solved in the usual way. Solved with 
respect to f 12 and with f 12 = A e w , a quadratic equation 
is obtained for oo 2 . With w - 1 \,- • 



2" C a x 2 + a/ ±7 (a, 2 - c^) + 4 a 2 2 a 4 2 ] 



(9) 



and we obtain for 



f 12 = A 1 sin X 1 t + A r cos \ x t + B x sin \ 2 t + B 2 cos \ 2 t 



i a 2 



H 2 + 



d t* 



12 
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* - An — sin \. t - Aa — — cos K- t 

X 2 — cc 3 X ^ 2 — cc 2 

~ Bt — r— — sin \, t - Bo — cos \ 2 t. 

The arbitrary constants are 



a 



Ap = 



Bi 



. ^L£?Ja. + (\ 2 _ a 3 ) d f 
2 a. t ^ 2 1 ; df 



( * j : i 

a 2 2 d t + (x i - a x } "TT^ 



x g ( \ 1 ° — \ 2 2 ) 



B 2 = -SIU. 



(X, 8 - \ 2 2 ) 



where 



d f i2 n , d f s3( 



1 12 0 > -23c' d~¥ " a " ia d t 

are the values for \ t -• 0. 

The initial conditions for t = 0 are 

d x 1 d Xj 

3E X = L. x + L 2 x 3 * 0 = - ca -g^t i; 0 

x 2 = l 2 _& - - c a 

d t 

when 

Xj_ - x 2 = L x - f 12 and x 2 - x 3 - L 2 - f 23 



12 o U ' *» 0 ~ U ' ' dt °» d t ~ ° a * 
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Hence 



A 2 = 0 3 2 - 0 



^ „ CC3 2 

"" 1 K ft?" - 57) ° a; 15 ^ " x 3 (x, 2 "- vi Ca 



For 

P a = kj i 12 and P 2 = kg f s3 
the forces on the hull or float are 

a 2 

p i = c a k i C 7 — I sin \ x t - 

\j ( \ ]_ - X 2 ) 



X 2 ( /v j 2 X 2 2 ) 



a 2 

sin X, t] (10) 



\ - °a *s [ ~ . V. . 1 , a , ^n X, t + 



X 2 - a 2 

+ fc sin \ 2 t] (11) 

A 2 \ A, j - A. s ; 



For a more convenient calculation^ let 

M = mass of the whole airplane - M a + M 2 , 

M x as mass of fuselage = rM, 

K 2 ss mass of float = sM, 

M 3 s= mass of water = wl£, 

k % = ck = elasticity between M x and M 2 , 

k 2 = ek s= elasticity between M g and M 3 , 



k - 



JL_JL- total elasticity between M x and M 3 . 



k x + k 2 
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Forces on the fuselage (mass ) : 

P x = G a /k M (qpi sin \ 1 t - cp 2 sin \ 2 t) (12) 

e c " eg 

s 



cp, - 7= ; cp 



By A + B 2 B,/A - 3 



1 [0 L+_5 + e UL5] 
<o r s s if 



X 



x- A (A+ 3); /| (A - 3) 



1 / / r + s s + w\ ,4ec 
.± / c - e • + x— 

2 J \ r s s w / s 2 



Forces on the float- (mass Mg): 

R = c aA /T: M (~ ^ sin \ x t + ty 2 sin \ 2 t) (13) 

f % a JL f A + B - C\ ^ c je_ /A - B - fl. \ 

3 B V A - B ? 8 2 3 V V ^A - B ; 



C = c 2L+JI 



9. When it is too difficult to calculate the spring con- 
stants, they can be determined by vibration tests on similar 
models. The seaplane is elastic ally suspended and caused to 
vibrate by a rotating weight eccentrically attached to the fuse- 
lage. Then, the natural vibration number of a structural member 
can be easily determined by its resonance with the revolution 
number of the eccentric weight , From this the spring constants 
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can be easily determined. A mean spring constant 
k = i/ki d Fi must be introduced for the elasticity of the 
float bottom, whence the mean bottom pressure is obtained. At 
the points where the spring constant of the bottom is smaller 
than the mean spring constant, the bottom pressure decreases, 
while it increases at the points where the constant is greater, 
namely, in the neighborhood of the bulkheads. 

10. A numerical example is given in the form of a calcula- 
tion for a float seaplane of the Heinkel monoplane type, e.g., 
the HE 5, HE 8, HE 9, etc. The elasticity between the fuse- 
lage and float is determined by a vibration test made by H. 
Hertel and Leiss of the Static Division of the D. V. L. on an 
HE 8 for the determination of wing vibration:.. 

The float was found to develop vibrations with a frequency 
of 550/min. , as shown in Figure 7. In the above equations 
therefore, we should put k^ = °°, since the axis of vibration 
passed through the rear suspension point of the float. However, 
as a first approximation, we shall use the formulas derived from 
the simplified assumptions which, strictly speaking, hold good 
only for k v = k h and a = b, or for a single spring located 
in the line of gravity of the constant k x • 

Let k x be the constant of a soring mounted in the line 
of gravity g h and producing the same type and frequency of 
vibration as that of the observed vibration. For the latter 



H.A.C.A. Technical Memorandum No* 580 24 

(Fig. 8) we have the approximations 

d 2 <p. 

CL 0 

(M s b ,2 + 0 2 ) = - k u e f y 

d t 

in which the following notation is used, 

M x , mass of fuselage = 275 kgs/m 

9 l5 inertia moment of fuselage = 825 kgm 2 

Mjg, mass of float = 30 kg 2 /m 

0 2 , inertia moment of float a 95 kgm 2 

Moreover, b = 0.8S and b' = 1.33. cp x and cp 2 are the 
angles of the fuselage and float motion; f v and k v the de- 
flection and spring constant of the forward suspension. For 

fv * e (c? x ~ cp 2 ) 



^ f U , ^ M r 1 



h& e 2 f p [ - + s ] 



d t* 1^ if + B t (Mg b 1 + G 2 ) 

We now replace the spring constant k v at a distance e by a. 

spring k x at a distance b and obtain the angular velocity 
of vibration 

</ 1 Vjii b 2 + 0j 



Mj b« 2 + 6 3 ^ 30 



By this formula kj = 574,000. 

This spring constant produces a deflection of 
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f 3700 kg a 7 mm 

f - 574000 kg/m - 4 *' 

under the static load of the weight of the fuselage. According 
to a numerical estimate, the spring constant at the bottom 
amounts to about k s /F » 1/750,000. For a float water line area 
of approximately 10 m 2 the mean deflection would be 

* 500 0 k S n <•«« 

f = — — "/ = 0.175 mm 

17500000 kg/m 

under the weight of the seaplane. 

The landing case represented in Figure 1 is used for the 
calculation. Let the airplane speed be V = 90 km/h = 35 m/s, 
the line of flight roughly horizontal and the angle of inclina- 
tion correspond to the angle of attack in leveling off (about 
12°). Hence, the normal speed component is o a = 5.2 m/s, the 
length of the bottom striking the water is approximately 1.3 m, 
according to Figure 1, and the spring constant of the bottom 
k 2 a 3,500,000 kg/m for a width of b = 0.36 m. 
For M = 305 kg 2 /m M x = 275 kg^/m 

S£ 2 =30 ' Mg = 2 | (a b 2 - |P) - 46 kg 2 /m 

fc 1 - 574000 § k 2 = 3500000 ^ 

k = * kg , = 483000 kg/m, 

(that is, r = 0.9; s * 0.1; w - 0.15; c - 1.16; e * 7.25), 

the impact on the fuselage is 

\ = 17800 sin 94 t - 3820 sin 450 t 
and the impact on the float is 
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p 2 = 568Q0 sin 450 t + 4820 sin 94 t. 

As shown by the calculation, the impact is a vibrational 
phenomenon. For the designing of airplanes it is therefore im- 
portant for the struts which, during the first moment of the 
impact, work in tension, to be subjected immediately to a com- 
pressive force of nearly the same magnitude. Some struts may 
also develop vibrations in resonance with the impact and there- 
fore collapse prematurely. The frequencies of the impact are 

n 1 = 15/s 
n s ~- 72/s* 

11. As w a s to be anticipated from the preliminary state- 
ments, the result was confined to the purely elastic impact 
with the impact coefficient 1, In practice the impact is 
damped, chiefly by the internal damping of the material and by 
friction in the connections, joints, etc. This damping action 
and its effect on the impact must be determined by tests. The 
question will be only briefly considered here. According to 
Plank? Honda and Koano, the damping of the material is directly 
proportional to the velocity of deformation, so that 

P = v f+ :pdf (15) 

If this formula were substituted for k f in equations 

(7) aid (8), the solution of the system fo r X s would gi ve a 

♦Plank, "3etrachtung uber dynamische Zugbeanspruchung , 11 Zeit- 
schrift des Vereines deutscher Irigenieure, 1912. 
**Honda and Konno, Zeitschrift fur arvrowandte Mathemat ik und 
Mechanik, 1921, p. 431. 
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biquadratic equation. According to the magnitude of the damping 
coefficient, the amplitudes thus obtained are smaller than those 
of undamped vibrations which, moreover, die out rather quickly. 
The damping coefficients must be determined experimentally. It 
is still uncertain whether the linear agreement is actually 
maintained. For plywood this assumption is not even approximate- 
ly correct. For metals the damping effect seems to reach values 
considerably above the proportionality limit, especially for 
duralumin. This is apparently the reason why the permissible 
load can be greatly exceeded with duralumin bottoms without 
causing failure, but merely bulging or other permanent deforma- 
tion. The above statement is based on the assumption that the 
permanent strength is not exceeded for the corresponding load- 
shifting coefficient and that the material is not impaired by 
corrosion. Owing to the brief duration of the process and to 
the damping effect above the limit of proportionality, the break- 
ing strength seems to be much greater than could be anticipated 
from the calculation based on static tensile tests. Within the 
elastic range, which alone is of interest here, the damping 
seems to be negligibly small, as shown by a short test made with 
a vibrating duralumin plate, and does not warrant the tedious 
calculation. It is therefore suggested, as the best approximate 
way of estimating the damping of duralumin floats, to use only 
the maximum value of the greater vibration for the maximum im- 
pact. The reason is that although, when damping is taken into 
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consideration, the amplitude does not differ much from that with 
no damping, the difference becomes apparent at a greater number 
of vibrations. Owing to the difference between the frequencies 
of the individual vibrations n 1 and n 2 , the maximum value of 
the impact force may be approximately expressed by the sum of 
the maximum values of the individual vibrations. With damping 
we can approximately assume instead that, when the maximum de- 
flection of the larger slow vibration is reached, the smaller 
but quicker vibration has already died out. Similar considera- 
tions apply to the impact force on a float. In this case the 
deflection of the short slow vibration remains small when the 
maximum value of the long, fast vibration is reached, while, for 
the maximum of the short slow vibration the long fast one has 
already largely died out. In this case the maximum value of 
the undamped vibration of great am;olitude can be substituted, 
with a fair approximation, for the maximum impact. Hence, the 
maximum impact force of the case calculated above is Pmiax^ 
17800 kg, or the load factor e « P/G = 5*9 and P 2 rnax = 
36800 kg or the bottom pressure 



~ "2 max _ — B 1#8 kg / cm 2 



F 2.06 

In spite of the limitation k v = these formulas can be 

used for the approximate calculation of a whole series of prob- 
lems, when mass and elasticity are properly subdivided* Thus, 
for float seaplanes, the impact force on the engine bearers or 
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wing can be calculated with a fair degree of accuracy, if the 
rest of the fuselage and the float are considered as a single 
mass Mg , if a mean elasticity between fuselage and water is 
put for k 2 and if the corresponding values of the engine bear- 
ers or wing are substituted for Ll 1 and k x . Similar calculations 
can be made for flying boats. A subdivision of the fuselage 
and wing is particularly advantageous for large flying boats, 
the weight of whose engines, fuel, etc., is distributed over 
the wings. 

12. In dividing a seaplane into two masses, the system of 
three combined vibrations mentioned in Section 7, which is rath- 
er difficult to calculate, is obtained for the eccentric impact. 
However, the problem seems to be covered sufficiently when the 
reduced mass is substituted for the mass M of the seaplane 
(equation 6) in the formulas for the step impact. On the sense 
of the calculation developed in this connection, the formula 
assumes the elasticity is located between the water and float 
and not distributed according to the above assumptions. 

13. The favorable influence of elasticity leads us to at- 
tempt a reduction of the impact by installing shock absorbers 
as on landing gears. This can be done to a certain extent. 

Yet too soft springs may easily have an effect contrary to the 
one desired. In taking off from a choppy sea and even from 
slightly rippling water, instabilities may develop under the 
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action of the accelerated mass of water subjected to short- 
period variations and owing to the small natural frequency of 
the system. Such instabilities, like the resonance of forced 
vibrations, may lead, to premature failure. This probably ac- 
counts for the failure of the repeated attempts to equip float 
gears with shock absorbers. 



14. Let us also investigate the V-shaped bottom for the 
case represented in Figure 1. A wedge of the length a (Fig. 
9), assumed to be without mass, is connected with the mass ¥i 1 
by a spring member. This wedge penetrates into the water at the 

d X 

time t with the speed ^ . A water mass Ivl, , no longer 
constant but a function of the width y, corresponds to the 
bottom portion of a width y and a depth x 2 immersed at the 
time t. The force on the bottom, which is assumed to have no 
mass, now equals the momentum increment of this variable mass 
with respect to time / r 



The V-Shaped Bottom 



P = 




d t 



As in the case of equation (3), we now have 




- k f 



> 



(16) 



d (u 



i d **\ 

d t 



= - k f 
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Moreover, % - x s = L - i and M 2 » f(y) » f(x 2 ). 

For very sharp V-bottoms and very large k values We can 

d x d x 2 d x . . . 

put — *= 3 ~ = -5 — r , and we then obtain 
1 a t a t d t ' 



H / v dxN 
d t ~ " d~t 



or 



d x _ y _ d X _j_ q ^ 



A,ai d t ~ " ^ & t 

The integration constant is 



dx 

unce in the case of t = 0 and M* =* 0. -rr = c„ 



0 = M 1 c a , 

" w > dt 



d x 1L 

= — ■ c n . 

d t Jl 1 + M a 

The impact force 

d 2 x li x s c a d M d x 



P 

As above, we again have 

d x 



[Mi + M 2 j d x d t 



1 « 

i-1 T M 2 



d t ~ Mi + ife 



a dy 

Moreover, tan — = -•— , vr hen the bottom walls are straight and 

2 dx 

make an angle of a with each other. Hence 



d It 



P - tan r - —or 



2 



2 ( li 1 + :i s ) 3 d y 

For the greatest width, which, with the usual bottom shapes, is 
likely to produce the maximum force, we have 
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P = tan * - a --3- 

2 / M^^ d y ( .. : _ b/2) 

V M a / 

With y - >,3 V _ 

Ma - P | (a b 2 - §-•) - | P (a y 2 - y 3 ), 

when a > "b 

d_H = ^ P f a h - | b»\ H 

d y . 3 V 4 / 

Y=h/3 

The impact force on the hull of a flying boat or on the 
fu3elage of a float seaplane with a sharp V-bottom is therefore 

P « tan § c a ? 1 b —4 (17) 

when a is the keel angle, 

c a the components of the landing speed normal to the keel, 

p the density of the water, 

a and b the length and width of the bottom striking the Water, 
a being smaller than b, 



;he mass of the seaplane, 

v / - 2 b 3N i -u 
^ - p f a d - ~- 1 when a > b. 

For twin-float or twin-hull seaplanes we would have 

tt pfa t> - -I b 8 ^ 

p . tttS o^ Si (is) 



ri + 



*A similar formula, though on the assumption of an infinitely 
long plate, was developed by Von Kariaan in a report published in 
October, 1929, by the National Advisory Committee for Aeronau- 
tics (Technical Note ho. 321: The Impact on Seaplane Floats 
during Landing). 
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in which 

M. = p f a o - 75-- I 

for a > To. 

Since, according to our assumptions, the sharp V- "bottom is 
not affected by elasticity, the development of a load factor ap- 
plicable to the whole airplane is warranted. The forces on the 
fuselage, as compared with its mass, are therefore smaller than 
the forces on the float bottom. 

According to the formulas, the impact is infinitely great 
for a = 180°. Even with a flat keel, the impact forces are 
excessive, so that the elasticity must then be taken into ac- 
count. The equation can be integrated numerically or graphically. 
As a first approximation, we might confine ourselves to deter- 
mining the impact of the flat bottom, taking elasticity into 
account, and the impact of the sharp V-bottom by the above formu- 
las. The impact of the flat-keeled bottom may then be approxi- 
mately determined by drawing through the point a = 180° a 
tangent to the curve of the sharp V- "bottom. 

We do not feel justified in expressing the keel of so-called 
wave-binding shapes (Fig* 10) by the angle a of Figure 10, as 
was hitherto usually done. In this case, equation (16) would 
have to be integrated graphically or numerically or else the im- 
pact of a flat bottom would have to be considered instead. 

15. The maximum bottom pressure is exerted on the immersion 
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of the keel. It is then 

p = tan c a 2 * p (19) 

for single-float or single-hull seaplanes, and 

p ;=£ tan | c a 2 n p (20) 

for twin-float seaplanes or twin-hull flying boats* The maximum 
stress on the bottom surface, however, is probably produced at 
a mean bottom pressure exerted over the whole bottom surface. 
This bottom pressure can be easily determined by formulas (17) 
and (18). 

IS. The eccentric impact can be calculated with the same 

i 2 

formulas if, as above, the reduced mass M' = M -r - is sub- 

i 3 + r & 

stituted for 1^ (page 15). 

Summary of the Results 

17. Figures 11 and 12 show the results which car* be theo- 
retically anticipated for Heinkel monoplanes calculated for a 
geometrically similar increase of dimensions but variations in 
landing speed and bottom angle. Inasmuch as the impact force 

for flat or Y-shaped bottoms is proportional to a surface, 

3 — jr 

P = oj G , can be expressed as a function of the weight. In 
this formula the coefficient c depends on the speed and the 
keel angle only. The load factor is then 
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The length of the waves which the seaplane can withstand 
is also increased, since the contact length of the bottom sur- 
face is included in the calculation and would also have to be 
proportionally increased, in order to preserve the geometrical 
similarity. 

The load factor calculated for the fuselage of the 3000- 
kilogram seaplane in seaway 2 was 6 g. It would be wrong, how- 
ever, to conclude that a seaplane calculated with this load fac- 
tor cannot resist stronger seaways, since a skilled pilot usu- 
ally succeeds in avoiding the case represented in Figure 1 by a 
tail landing. Yet seaway 2 seems to be the limit at which a 
pilot can bring his plane down without special training. Be- 
sides, seaways seldom correspond to conditions which can be rep- 
resented diagrammat ically. -However, such a representation is 
also used in shipbuilding practice for strength calculations 
and is necessary in order to obtain a basis for the calculation. 

Comparison with the DVL Load Assumptions 

18. The load, assumptions of the D f V.L« developed from data 
supplied by Lewe and the experience of various companies are 
based on geometrical similarity. The load assumptions do not 
account for the influence of bottom width, mass distribution, 
elasticity, etc* Figure 13 shows the influence of the flying 
weight on the load factor, according to theory and to the D.V.L. 
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load assumptions.* Figure 14 illustrates the influence of the 
bottom angle on the impact according to the two methods of cal- 
culation* There is also plotted a point taken from Bottomley's 
model tests,** which is closer to the theoretical values than 
to those of the D«V.L« load assumptions. According to theory, 
the velocity is proportional to the first power for flat bottoms 
and to the second power for sharp V-bottoms, while it is propor- 
tional to the 1.5 power under the D.V.L. load assumptions. 

According to these load assumptions the bottom pressure is 
calculated from the 50;£ greater fuselage loading. In theory 
the float has about twice the fuselage loading, provided it is 
not carrying additional loads (fuel tanks). The area over which 
the load is distributed is about the same in theory and accord- 
ing to the load assumptions. In general, it can be said that 
the theory is not in fundamental contradiction with the empir- 
ically developed load assumptions. One advantage of the theory 
over the load assumptions lies, however, in the possibility of 
considerably more accurate calculations and thus hitting the 
best compromise of the different float factors, especially as 
regards quick take-off and adequate strength. The theory, how- 
ever, requires experimental confirmation and extension by experi- 

ments which are now under way and which w ill soon be reported. 

•Lewe, Zeitschrift fur Flugtechnik und Motorluf tschif f ahrt , 
1920, p. 125. 

**Bottomley, "The Impact of a Model Seaplane Float on Water." 
British A.G.A. Reports and Memoranda No. 583 (1913). 
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Fig.l Seaplane alighting in seaway 2. 
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Figs. 1,4, 5 
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Fig. 4 Accelerated '.water mass with one sided flow 



Fig, 5 Diagram for the seaplane float system, 
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Figs. 2 ,3 
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Fig. 6 



Diagram of eccentric impact cf the two-mass 
system . 
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Fig. 7 Vibration test for the determination of the 

elasticity between fuselage and float , carried 
out on a Heinkel HE 8 seaplane. 
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Fig. 8 Diagram of the test according to Fig. 7, 
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Figs. 9, 10,11 
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Fig. 9 Diagram for V bottom tests. 




Fig. 10 V bottom of vrav e -binding 
shape . 
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Fig. 11 Alighting impact forces of twin-float 
seaplanes of the Heinkel monoplane 
type (determined theoretically). 
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Figs. 12, 13,14 



a, V=130 km/h 

b, V=110 km/h 
o,V= 90 km/h 
d,v= 70 km/h 
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Fig. 12 Water pressures at the step of t^rin-float seaplanes 
of the Heinkel monoplane type. (p n = f [oc f V] 

theoretically) . The pressures are independent of the 
increases In size of the seaplanes. 
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Fig. 13 Relation between load factor and flying- weight 
of similar aircraft. 
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